EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 



CERN-EP/2000-157 

1st March 2000 



Measurement of the semileptonic b 
branching fractions and average b 
mixing parameter in Z decays 



DELPHI Collaboration 



Abstract 



The semileptonic branching fractions for primary and cascade b decays 
BR(6^£^), BR(b— >c^£+) and BR(b^c^£^) were measured in hadronic Z 
decays collected by the DELPHI experiment at LEP. 

The sample was enriched in b decays using the lifetime information and various 
techniques were used to separate leptons from direct or cascade b decays. 
By fitting the momentum spectra of di-leptons in opposite jets, the average b 
mixing parameter x ^ilso extracted. 
The following results have been obtained: 



BR(6-^r) = {10.70 ±0.08{stat)±0.2l{syst)-°ilimodel))% 

BR(b^c-^£+) = {7.98±0.22{stat)±0.21{syst)to,ltimodel))% 

BR(b^c^r) = {1.6l±0.20{stat)±0.17{syst)tofi{model))% 

X = 0.127 ± 0.013(stot) ±0.005(syst) ±0.004(moc/e/) 



(Accepted by Eur.Phys.J.C ) 
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1 Introduction 

Measurements of the direct semileptonic branching fractions of 6-hadrons are impor- 
tant in order to understand the dynamics of heavy quark decays and to determine the 
weak couplings of quarks to the W boson. From a precise measurement of the inclusive 
semileptonic branching fractions of b quarks a precise value of the Cabibbo-Kobayaski- 
Maskawa matrix element \Vcb\ can be calculated 

These measurements have been performed both at the T(4S') and in hadronic Z de- 
cays. In order to make a comparison between the two sets of results, the fact that the 
composition of the inclusive sample is different in the two cases must be taken into ac- 
count. At low energy only B~ and mesons are produced, while at the Z , mesons 
and 6-baryons are also present. Assuming the semileptonic widths of all 6-hadrons to be 
equal, their respective semileptonic branching fractions are expected to be proportional 
to their measured lifetimes. The ratio between the B~ and B'^ lifetimes to the inclusive 
6-hadron lifetime measured at the Z , is at present larger than 1, whereas the semilep- 
tonic branching fractions of 6-hadrons measured at the Z are slightly larger than the ones 
measured at the T(45) . 

Theoretical calculations which include higher order perturbative QCD corrections give 
a prediction of the branching fraction value correlated with the prediction for < Uc >, 
the average number of charmed hadrons produced per 6-hadron decay These results 
are compatible with the present LEP measurements. 

In this paper, the two cascade processes: 6— >c— and b^c-^i~ are also consid- 
ered, not only because they are the main source of background to the direct decays, 
but also because the values of these branching fractions are important inputs to several 
other heavy flavour measurements, like asymmetries and oscillations measurements. The 
BR(b-^c-^^~) measurement presented in this paper is the first inclusive measurement of 
"right sign" leptons from cascade decays of 6-hadrons. 

In addition, the average B^ — mixing parameter is measured. It is the time in- 
tegrated probability that a 6-hadron oscillates into a 6-hadron: x = ^^^"f^Zlfix^'^^ ■ It 
is related to the mixing parameters of 5° and B^ mesons, Xd and Xs respectively, by: 
X = 9B°Xd + 9boXs, where g^o and gso are the production fractions of B^ and S° in 
semileptonic decays. Its measurement can therefore be used in the evaluation of the 
production fraction of 5° mesons ||l|. 

This paper presents the measurement of inclusive semileptonic branching fractions of b 
quarks in hadronic Z decays using data collected with the DELPHI detector at LEP. Four 
analyses have been performed, using different strategies and using various data samples, 
partially overlapping. Events containing b hadrons were selected using lifetime informa- 
tion, electrons and muons were identified and several different techniques were used to 
determine the origin of the lepton. Direct and cascade branching fractions: BR(6— 
BR(b—>c— and BR(b^c^£~) were measured and, by fitting the momentum spectra 
of di-leptons in opposite jets, the average B^ — B^ mixing parameter x was also extracted. 

The previous DELPHI results on the semileptonic branching fractions ||^ were ob- 
tained with data collected at LEP in 1991 and 1992, using electrons and muons in a 
sample of hadronic Z decays, with natural composition of quark fiavours. A global fit 
to several electroweak parameters was performed. With respect to that analysis there 
is little dependence on the partial decay widths of the Z into bb and cc quark pairs 
{Rb = Tjji/Thad, Rc = ^ccl^had) and the background due to misidentified hadrons and 
leptons from decays and punch-through of light hadrons has been reduced. The present 
result supersedes the previous result obtained by DELPHI ||5|. 
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The layout of the paper is the following: a description of the DELPHI detector is given 
in Section ||. The selection of the hadronic event sample is described in Section ^ The b- 
flavour tagging algorithm is described in Section ^. A brief summary of the performances 
of lepton identification algorithms is given in Section |. Results obtained in the different 
analyses are then described in the following Sections: single and di-lepton analysis (Sec- 
tion H), single lepton and jet-charge analysis (Section |^), multitag analysis (Section ||) 
and inclusive 6-hadron reconstruction analysis (Section D . Finally, in Section |1^ averages 
of the results obtained in the different analyses are calculated. 

2 The DELPHI detector 

The DELPHI detector has been described in detail in reference P|. Only the compo- 
nents relevant to this analysis are mentioned here. 

In the barrel region, the charged particles are measured by a set of cylindrical tracking 
detectors with a common axis parallel to the 1.2 T solenoidal magnetic field and to the 
beam direction. The time projection chamber (TPC) is the main tracking device. The 
TPC is a cylinder with a length of 3 m, an inner radius of 30 cm and an outer radius of 122 
cm. Tracks are reconstructed using up to 16 space points in the region 39° < 9 < 141°, 
where 9 is the polar angle with respect to the beam direction. Tracks can be reconstructed 
using at least 4 space points down to 21° and 159°. 

Additional precise i?$ measurements, in the plane perpendicular to the magnetic field, 
are provided at larger and smaller radii by the Outer and Inner detectors, respectively. 
The Outer Detector (OD) has five layers of drift cells at radii between 198 and 206 cm 
and covers polar angles from 42° to 138°. The Inner Detector (ID) is a cylindrical drift 
chamber having inner radius of 12 cm and outer radius of 28 cm and covers polar angles 
from 23° to 157°. It contains a jet chamber section providing 24 i?$ coordinates measure- 
ments surrounded by five layers of proportional chambers with both and longitudinal 
z coordinates measurements. 

The micro-vertex detector (VD) is located between the LEP beam pipe and the 
ID. It consists of three concentric layers of silicon micro- vertex detectors placed at radii 
of 6.3, 9.0 and 10.9 cm from the interaction region, called closer, inner and outer layer, 
respectively. For all layers the micro-vertex detectors provide hits in the i?<I>-plane with 
a measured intrinsic resolution of about 8 /im; the inner and outer layers provide in 
addition measurements in the z direction, with a precision depending on the polar angle 
and reaching a value of 9 /xm for tracks perpendicular to the modules. The polar angle 
coverage for charged particles hitting all three layers of the detector is 44° < 9 < 136°; 
the closer layer coverage goes down to 25°. The z measurement was only available in 
1994 and 1995. 

Additional information for particle identification is provided by the Ring Imaging 
Cherenkov counters (RICH) measuring the Cherenkov light emitted by particles travers- 
ing a dielectric medium faster than the speed of light. The barrel part of the detector 
covers polar angles from 40° to 140°. To cover a large momentum range, a liquid (C6-F14) 
and a gas (C5F12) radiator are used. 

The barrel electromagnetic calorimeter, HPC, covers the polar angles between 42° 
and 138°. It is a gas-sampling device which provides complete three dimensional charge 
information in the same way as a time projection chamber. Each shower is sampled nine 
times in its longitudinal development. Along the drift direction, parallel to the DELPHI 
magnetic field, the shower is sampled every 3.5 mm ; in the plane perpendicular to the 
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drift the charge is collected by cathode pads of variable size, ranging from 2.3 cm in the 
inner part of the detector to 7 cm in the outer layers. 

In the forward regions the tracking is completed by two sets of planar drift chambers 
(FCA and FCB) placed at distances of ±165 cm and ±275 cm from the interaction point. 
A lead glass calorimeter (EMF) is used to reconstruct electromagnetic energy in the 
forward region. 

For the identification of hadronic showers, the iron return yoke of the magnet is in- 
strumented with limited streamer mode detectors to create a sampling gas calorimeter, 
the Hadronic Calorimeter (HAC). 

Muon identification in the barrel region is based on a set of muon chambers (MUB), 
covering polar angles between 53° and 127°. It consists of six active planes of drift 
chambers, two inside the return yoke of the magnet after 90 cm of iron (inner layer) and 
four outside after a further 20 cm of iron (outer and peripheral layers) . The inner and 
outer modules have similar azimuthal coverage. The gaps in azimuth between adjacent 
modules are covered by the peripheral modules. Therefore a muon traverses typically 
either two inner layer chambers and two outer layer chambers, or just two peripheral 
layer chambers. Each chamber measures the i?$ coordinate with a precision of about 2- 
3 mm. Measuring /?$ in both the inner layer and the outer or peripheral layer determines 
the azimuthal angle of muon candidates leaving the return yoke within about ±1°. These 
errors are much smaller than the effects of multiple scattering on muons traversing the 
iron. 

In the forward region the muon identification is done using two sets of planar drift 
chambers (MUF) covering the angular region between 11° and 45°. The first set is placed 
behind 85 cm of iron and the second one behind an additional 20 cm. Each set consists 
of two orthogonal layers of drift chambers where the anode is read out directly and the 
cathode via a delay line to measure the coordinate along the wire. The resolution in both 
coordinates is about 4 mm. 

3 Event selection 

Charged particles were accepted if their polar angle was between 20° and 160°, their 
track length was larger than 30 cm, their impact parameter relative to the interaction 
point was less than 5 cm in the plane perpendicular to the beam direction and less than 
10 cm along the beam direction and their momentum was larger than 200 MeV/c with a 
relative error smaller than 100%. Neutral particles detected in the HPC and EMF or in 
the hadronic calorimeters were required to have a measured energy larger than 500 MeV. 

The decays of the Z to hadrons were selected by requiring a total energy of the charged 
particles (assumed to be pions) larger than 15% of the center-of-mass energy and at least 
7 reconstructed charged particles. With these criteria, the efficiency to select qq events 
from the simulation was about 95%. All sources of background have been found to be 
below 0.1%. No significant differences in the acceptance between different fiavours have 
been found. 

For each event the thrust axis was calculated from the selected charged and neutral 
particles. Only events with: \co^dthrust\ < 0.90 were used. Requiring, in addition, that 
all sub-detectors needed for these analyses were fully operating, totals of about 1 030 000 
and 515 000 Z hadronic decays were selected from the 1994 and 1995 data samples, 
respectively. About 3 800 000 events were selected from a simulated sample of Z ^ qq 
events. A reduced angular region was used in some parts of the following analyses to 
ensure an efficient acceptance for the vertex detector. 
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Events were generated with the JETSET 7.3 generator [|] using parton shower and 
string fragmentation with parameters optimized to describe the hadronic distributions as 
measured by DELPHI 0. Generated events were passed through a detailed simulation 

which modeled the detector response and processed through the same analysis chain 
as the real data . Jets were formed from the charged and neutral particles using the 
JADE algorithm with Y^l^ = 0.02 ||T0|. The transverse momentum of the lepton ( pt ) 
was determined relative to the direction of the jet, excluding the lepton itself. 

Any differences with respect to these selection criteria, as well as their effect on the 
statistics used, will be explicitly described for each analysis. The four analyses used 
different data subsamples corresponding to the optimal operation of the subdetectors 
relevant to the definition of the variables used. Analysis I and IV used 1994 and 1995 
data samples, Analysis III used also 1992 and 1993 data, while Analysis II used 1994 
only. The 1992 and 1993 statistics are given in Section ||. 

4 b-flavour tagging 

A 6-flavour tagging algorithm was used in order to obtain a sample enriched in Z ^ bb 
events. Events were divided into two hemispheres, with respect to a plane perpendicular 
to the thrust axis and passing through the beam interaction point. The fe-flavour tagging 
algorithm was applied separately to each hemisphere. Analyses I and IV used the com- 
bined 6-flavour tagging algorithm described in ||ll|. This algorithm combines, in a single 
variable, several quantities which are sensitive to the presence of a 6-hadron. 

The main discriminant variable is the probability for all tracks belonging to the hemi- 
sphere to come from the primary vertex, calculated from the impact parameters of the 
tracks positively signed according to the lifetime convention. Other variables were defined 
for hemispheres containing a secondary vertex. These variables are: the effective mass 
of the system of particles attached to the secondary vertex, the rapidity of these tracks 
with respect to the jet direction and the fraction of the charged energy of the jet which 
is included in the secondary vertex. Optimized levels of efficiency and purity were chosen 
in each analysis. 

Analysis II used a 6-fiavour tagging algorithm exploiting only the information from 
the impact parameters of charged particles ||ll|. Analysis III used a multivariate method 
to tag the fiavours, as described in Section ^T] . 

5 Lepton sample 

5.1 Muon identification 

To identify a charged particle with momentum greater than 3 GeV/c as a muon can- 
didate, its track was extrapolated to each of the layers of the muon chambers taking into 
account multiple scattering in the material and the propagation of track reconstruction 
errors. A fit was then made between the track extrapolation and the position and direc- 
tion of the hits in the muon chambers. Ambiguities with muon chamber hits associated 
to more than one extrapolated track were resolved by selecting the track with the best 
fit. The charged particle was then identified as a muon if the fit was sufficiently good and 
if hits were found outside the return iron yoke. 
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To exclude regions with poor geometrical acceptance, a muon was accepted only if its 
polar angle, 6^, was within one of the following intervals: 

0.03 < I cos^^l < 0.62 or 0.68 < | cos^^| < 0.95, 

which defined the barrel and the forward regions, respectively. 

The muon identification efficiency was measured in Z — > events, in the decays 

of taus into muons and using muons from two-photon collisions 77 — > yU'*'/^". A mean 
efficiency of 0.82 ± 0.01 was found with little dependence on the muon momentum and 
on the track polar angle. Predictions of the simulation agree with corresponding mea- 
surements in data, both in absolute value and in the momentum dependence, within a 
precision of 1.5%. 

An estimate of the misidentification probability was obtained by means of a lifetime- 
based anti 6-tag to select a background enriched sample. After the subtraction of the 
muon content in the selected sample, the misidentification probability was found to be 
(0.52 ± 0.03)% in the barrel and (0.36 ± 0.06)% in the forward regions. Applying the 
same procedure to the simulation gave however lower values, with factors 2.03 ± 0.12 
(2.02 ± 0.13) in the barrel and 1.22 ± 0.20 (1.78 ± 0.24) in the forward regions for the 
1994 (1995) samples, respectively, showing a small momentum dependence and about 
30% reduction near the borders of the geometrical acceptance of the muon chambers. 

The hadron misidentification probability, measured both in data and in simulation, 
was cross-checked using pions from and r decays and compatible results were found. 
In Analysis I, II and IV the simulated hadrons misidentified as muons were reweighted 
according to the probability measured in data. In Analysis III a different approach was 
used to estimate the misidentification probability, as described in Section |87^ , and good 
agreement with the above results was found. 

5.2 Electron identification 

Charged particles with momenta greater than 3 GeV/c and within the efficient accep- 
tance region of the HPC (0.03 < | cos6'e| < 0.72) were selected as electron candidates on 
the basis of the information from the HPC, the TPC and the RICH detectors. Tracks 
were extrapolated to the HPC and associated to detected showers. The signals from the 
various detectors were then analyzed by a neural network. By using the network response 
obtained in a sample of simulated electrons from b and c decays, a momentum dependent 
cut was defined in order to have a 65% efficiency, constant over the full momentum range. 

To reduce the contamination from electrons produced from photon conversions, elec- 
tron candidates were removed if they came from a secondary vertex and carried no trans- 
verse momentum relative to the direction from the primary to this secondary vertex. 

The efficiency of tagging an electron was measured in the data by means of a sample 
of isolated electrons extracted from selected Compton events and a sample of electrons 
produced from photon conversions in the detector. The ratio between the values of the 
efficiencies measured in real and simulated events was parameterized in terms of the pt 
and the polar angle of the track and found to be on average 0.92 ± 0.02 and 0.93 ± 0.02, 
in the 1994 and 1995 samples, respectively. A corresponding correction factor was then 
applied to the sample of electrons in simulated qq events. 

The probability of tagging a hadron as an electron was also measured in the data by 
selecting a background sample by means of the anti 6-tag technique in the same manner 
as for muons. The measured misidentification probability in data and the ratio with the 
same quantity obtained in simulated events were on average (0.40±0.02)% and 0.76±0.05 
in the 1994 sample and (0.38 ± 0.04)% and 0.70 ± 0.06 in the 1995 sample. 
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5.3 Simulated lepton sample 

Samples of simulated events, which were processed through the same analysis chain as 
the data as described in Section ||, were used to obtain reference spectra for the different 
sources of simulated leptons. 

The b semileptonic decays to electrons and muons were simulated using the model of 
Isgur et al. ||l2j (ISGW model in the following). The model of Bauer et al. which 
takes into account the finite mass of the produced lepton, was used for the b decays into 
r's. For D decays the branching ratios were adjusted to be in better agreement with 
measured values f^. In the different semileptonic decay modes, the branching fractions 
for the decays to neutral pions, when not measured, were obtained imposing isospin 
invariance. Reference spectra with alternative models have been obtained reweighting 
the events according to the decay model considered. The weight was computed on the 
basis of the lepton momentum in the B{D) rest frame. According to the prescription 



of IT^, for the central value of the results, the inclusive model of Altarelli et al. ||T5 
(ACCMM model in the following) was used, with model parameters tuned to the CLEO 
data ||16|, whereas ISGW and ISGW** models have been used to evaluate the systematic 
uncertainties. ISGW** indicates the ISGW model modified to include a 32% contribution 
of charmed excited states (referred to as D**), instead of the original 11% predicted by 
the model itself, so as to better describe the CLEO data. 

Leptons from the decay chain b — ^ cW — > ccq — ^ ci~X (the so called "upper decay 
vertex") were considered with the contributions from both Dg i~X and D^{D~) 

i-x. 



6 Analysis I: Measurement of semileptonic b decays 
from single leptons and di-leptons spectra 

In this analysis the semileptonic branching fractions for primary and cascade b decays 
BR(6^£~), BR(b^c^£+) , BR(b^c^£~) and the average b mixing parameter, are 
measured using the momentum spectra of single lepton and di-leptons in opposite jets. 
The single lepton spectra are studied in a sample of events highly enriched in bb, selected 
by means of a 6-fiavour tagging algorithm. In the di-lepton sample, the bb purity is 
increased by requiring a minimum pt for one of the leptons. 

The sensitivity to the different sources of leptons is given by the kinematic properties 
of leptons from different sources and by the charge correlation between di-leptons in 
opposite jets from b and b, respectively. 

Hadronic events and lepton candidates were selected as described in Sections || and |5|. 
The angular region | cos 9thrust\ < 0.9 was used for di-lepton candidates, while for single 
lepton events, to have a good efficiency in the fe-flavour tagging, events were considered 
only if they fulfilled \cos9thrust\ < 0.7. As a consequence, only barrel muon chambers 
were considered for single muons. About 768 000 and 385 000 Z hadronic decays were 
selected in the 1994 and 1995 data samples, respectively. 



6.1 Single lepton fit 

Events were divided into two hemispheres with respect to a plane perpendicular to the 
thrust axis and passing through the beam interaction point. A primary vertex was recon- 
structed in each hemisphere to suppress possible correlations between the two hemispheres 
induced by the 6-tagging algorithm. The combined 6-fiavour tagging algorithm described 
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in Section |^ was used to select hemispheres enriched in 6-hadron content while, in the 
opposite hemisphere, the single lepton spectra were studied. For the cut on the combined 
6-tagging variable used in this analysis, the following efficiencies for selecting different 
flavours were estimated from simulation: £5 = (39.34 ± 0.05)%, Sc = (1.87 ± 0.02)%, 
^uds = (0.189 ± 0.003) %, so that the fraction of b events in the sample was Vb = 95.1%. 

The value of Eb is quoted only for reference, since it is never used in the following. In 
practice the number A^"^ of tagged hemispheres which contain a b quark was estimated 
as: 

= Nt^ag - i^c y<Rc + Suds X Ruds) X 2Nhad 

where: N^^ and Nhad are the total numbers of tagged hemispheres and the number of 
hadronic events, respectively, Ec and Euds were the efficiencies for charm and light quark 
events, respectively, obtained from simulation, and Ruds = ^uds/^had = ^ — Rb — Rc- 
The LEP averages of 0.21643 ± 0.00073 and 0.1694 ± 0.0038 were used for Rb and R^, 
respectively The number of bb events used in the simulation was normalized to the 

same value Nf^ 

Once a hemisphere was tagged as b, leptons were studied in the opposite hemisphere. 
A correction was applied, estimated from simulation, because of the correlation between 
the lifetime and the lepton tags. It arose mainly from the acceptance requirements, 
which are different for electrons and muons, and amounted to pe = 1.003 ± 0.005 and 
Pij, = 1.017 ± 0.005. Here p is the fraction of lepton candidates found in the hemisphere 
opposite to the 6-flavour tagged hemisphere, compared to the fraction of lepton candidates 
found in an unbiased b hemisphere. Before calculating the lepton transverse momentum, 
a search for secondary vertices was performed using the same algorithm as in ||1 1|| . 
When the secondary vertex was successfully reconstructed (about 45% of the events), the 
primary to secondary vertex direction was found to give a better approximation of the 
6-hadron flight direction than the jet axis, and was used in its place. The resolution on 
the 6-hadron flight direction improved correspondingly from 30 to 20 mrad. 

Lepton candidates were classified according to their different origin as follows: 

direct 6-decay: 

b^i- + X, 

"right sign" cascade decays: 

b^c + x^e- + x, 

"wrong sign" cascade decays: 
b^c + X^i-^ + X, 
b decays into r lepton: 

b^T- + x^e- + x, 

direct c-decay 
c ^ + X, 

prompt leptons from J/\i/ decays or from 6 or c decays, where the cc (bb) pair is 
produced by gluon splitting, 
misidentified or decaying hadrons. 

The above classification was considered both for electrons and muons, separately. 

A binned maximum likelihood fit was used to compare the momentum and transverse 
momentum spectra of electrons and muons in data with the simulation. The full likelihood 
expression is reported in appendix. 
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6.2 Di-lepton fit 

The single lepton likelihood was multiplied by a likelihood obtained for di-leptons in 
opposite hemispheres, in order to separate the b^i~ from the 6— >c— and the b^c^£~ 
components and to extract the average mixing parameter x- In the di-lepton sample no 
6-flavour tag was used in order not to introduce any bias in the composition of the b- 
hadron sample. The b enrichment was obtained by requiring a minimum pt for one of the 
two leptons. The full pt spectrum was considered for the opposite lepton. For a cut at 
Pt > 1.2 GeV/c, a b purity of about 88% was obtained using simulated events. 

Di-lepton events were separated, for both the data and the simulated samples, into six 
groups depending on whether the two lepton candidates have the same or opposite charge 
and on which combination of lepton species (ee, efi, ii/i) they belonged to. Lepton pairs 
were used if the two leptons were separated by at least 90", while lepton pairs coming 
from the same jet were omitted from the fit to avoid additional systematic uncertainties 
in the composition of the cascade lepton sample. In each group, simulated events were 
separated into di-lepton classes, according to the different possible combinations in the 
two hemispheres of the above mentioned single-lepton classes (a) to (g). To guarantee 
a reasonable number of events in each bin, the p and pt of each lepton in the pair were 



combined to form a single variable, the combined momentum, Pc, defined as in 19 



Pc = y Pt ^ ■ Two-dimensional reference distributions were obtained for the chosen 

combinations in the variables where (p"^""^) refers to the smaller (larger) 

combined momentum. 

If — mixing is not considered, the main source of di-leptons having opposite 
charges are direct 6-decays: (5 — ^ £~)(p ^''')- But, in the presence of mixing, a fraction 
2x(l~x) of these di-leptons have the same charge. Same charge di-leptons also arise from 
events with one direct 6-decay and one cascade 6-decay: (6 — > ^ c ^ Because 
of mixing, a fraction 2^(1 — x) of these events will enter the opposite charge class. 

The fraction of leptons of class a, b and c were determined by the fit, whereas con- 
tributions from lepton classes (d) to (g) were fixed to the values given in Table |^. The 
detailed expression of the likelihood function, for single lepton and di-lepton, is reported 
in appendix. 



6.3 Results and systematic uncertainties 

The results obtained with the 1994 and 1995 samples and their average are shown in 
Table [l|, where the uncertainties are statistical only. About 12% of the single leptons were 
also included in the di-lepton sample and the statistical uncertainties have been corrected 
accordingly. 





1994 


1995 


1994±1995 


BR(6^r) 
BR(b^c->£+) 
BR(b^c^r) 
X 


0.1066 ±0.0014 
0.0822 ±0.0049 
0.0144 ±0.0044 
0.119 ±0.016 


0.1081 ±0.0019 
0.0781 ±0.0064 
0.0196 ±0.0056 
0.138 ±0.022 


0.1071 ±0.0011 
0.0805 ±0.0039 
0.0164 ±0.0035 
0.126 ±0.013 



Table 1: Results of the fit to the 1994 and 1995 lepton samples and their combination. 
The uncertainties are statistical only. 
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In Figure |l] single lepton and di-lepton spectra are shown. The simulation spectra 
have been reweighted according to the result of the fit. The correlation matrix for the 
statistical uncertainties is shown in Table 0. 





BR(6^r) BR(b^c^£+) BR(b^c^r) x 


BR(6-^r) 
BR(b^c^^+) 
BR(b^c-^r) 
X 


1.00 -0.241 -0.061 0.086 
1.00 -0.797 -0.159 
1.00 0.112 
1.00 



Table 2: Correlation matrix of statistical uncertainties in Analysis I. 

The following sources of systematic uncertainties have been considered: 

• experimental uncertainty related to lepton measurements: 
the muon and electron identification efficiencies and the background due to hadron 
misidentification have been varied considering their measurement uncertainties in the 



data-simulation comparisons (see Sections |5T1|J5.2| ) . To account for effects related to 



the difference in topology between the test samples used in Sections |5Tl| , |5^ and 
the hadronic environment, an additional uncertainty of ± 2% has been applied to 
the efficiencies, as estimated from simulation. As a consequence, the total relative 
uncertainties assumed on the leptons efficiencies were ± 2.5% and ± 3% for muons 
and electrons, respectively. The residual contamination in the electron sample due 
to converted photons has been varied by ± 10%. 

The angular distribution between di-leptons is well described by simulation, therefore 
the angular cut of 90° is assumed not to add any systematic uncertainty. 
The fit has been performed using for the pt calculation both the jet direction and the 
secondary vertex direction. Half the difference between the results has been used as 
systematic uncertainty. 

experimental uncertainty related to the 6-fiavour tagging: 

efficiencies to tag c and uds quarks have been varied by 9% and 22%, respectively, 
according to the uncertainties in ||ll|. The partial decay widths Rh and Rc have 
been varied according to their measurement uncertainties. 

The correction factors for the correlation between the 6-tag and the leptons (pe , Ppi) 
have been varied by twice their statistical uncertainties. The dependence on lepton 
momentum of the correlation has also been studied. Since the 6-tag efficiency is 
higher in presence of high momentum leptons, the lepton spectrum in hemispheres 
opposite to a 6-tagged one is slightly biased towards low momenta. A correction has 
been estimated with simulation comparing spectra in tagged and non tagged events 
and the full effect has been assumed as a systematic uncertainty. 
The stability of the result as a function of the cut on the 6-fiavour tagging variable 
has been checked to be compatible with the corresponding statistical fiuctuations. 
modelling uncertainty related to the assumed physical parameters: 
the mean value and the range of variation of several physical parameters used in 
the simulation was calculated according to references [|], ||l4l and [|l3]. In particular 



they have been varied: the mean fractional energy of h and c hadrons, the branching 
fractions assumed for 6 ^ r — ^ £, 6 — > J/\E' — > c — *• £ and the fraction of gluon 
splitting to heavy quarks. The lepton distribution from the "upper vertex" was 
studied by varying the contributions of Ds i^X and D^[D^) i^X of the 
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Figure 1: Comparison of data and simulation spectra. The simulation spectra have 
been reweighted according to the result of the fit. (a) Transverse momentum distribution 
for single electrons and muons. b — > x indicates h decays to misidentified or decaying 
hadrons. (b)((c)) Combined momentum distribution for the two leptons in di-lepton 
events, identified in opposite jets and having the opposite (same) charge, p™" refers to 
the minimum combined momentum of the two leptons. In the legend of (b) and (c) the 
lepton origin in the two hemispheres is described, the label "mix" refers to events where 
mixine occurred. 



Error Source 


Range 


ABR(6^£-) 
10-2 


ABR(b^c^£+) 
10-2 


ABR(b^c^r) 

10-2 


Ax 

10-2 


electron efficiency 




±3% 


^0.15 


^0.14 


TO.06 


±0.02 


misidentffied e 




±8% 


TO.05 


TO. 14 


TO.06 


±0.04 


converted photons 




±10% 


<0.01 


TOM 


TO.03 


±0.01 


/i efficiency 




±2.5% 


to.u 


TO. 18 


TO.05 


±0.06 


misid. /i barrel, forward 


±6.5%, 17% 


TOM 


TO. 15 


TO.06 


±0.02 


jet direction 




see text 


±0.05 


-0.03 


-0.08 


+ 0.6 


£c 




±9% 


±0.02 


TO.01 


TO.01 


±0.03 


^uds 




±22% 


±0.01 


±0.02 


<0.01 


TO.02 


i — b correlation 




±1% 


^0.05 


TO. 11 


TO.03 


±0.03 


i — b corr. p dependence 


see text 


=F 0.04 


± 0.03 


TO.01 


T 0.04 


Rb 




0.21643 ± 0.00073 [17] 


<0.01 


<0.01 


<0.01 


<0.01 


Rc 




0.1694 ± 0.0038 [17] 


<0.01 


<0.01 


<0.01 


<0.01 


XE{b) 




0.702 ± 0.008 [14] 


=F0.11 


±0.07 


±0.04 


TO. 15 


xe{c) 

b^W^D 
b^W^Ds 

BR(6 -^T-^t) 




0.484 ± 0.008 [14] 

(l-28lo:6l) [14] 
(0.459 ±0.071)% [2] 


=F0.02 
±0.03 
^0.02 


±0.03 

+0.20 
-0.11 

TO.03 


TO.03 

—0.23 
+0.13 

TO.04 


±0.02 

—0.09 
+0.07 

±0.02 


BR(6 ^ J/^Ij t 


n 


(0.07 ± 0.01)% [2] 


=F0.03 


±0.01 


±0.01 


TO.09 


BR(c^£+) 




f Ci o _i_ n QO^0/ Fi ^1 
[y.oD ± U.ozj/o [i l\ 


=FO.Oi 


TO. 03 


TO.04 


1 r\ n 1 

±0.0i 


g ^ cc 




(3.19 ±0.46)% [17] 


<0.01 


<0.01 


<0.01 


<0.01 


g ^ bb 




(0.251 ± 0.063)% [17] 


<0.01 


<0.01 


<0.01 


±0.01 


total systematic 




±0.26 


±0.38 


±0.25 


±0.64 


Semilept.mod.b 
Semilept.mod.c — » 


414] 


ACCMM (+S„) 
ACCMM1(11™) 


-0.24 
+0.41 
-0.08 
+0.07 


+0.23 
-0.29 
-0.11 
+0.01 


+0.14 
-0.23 
-0.03 
+0.02 


-0.23 
+0.28 
-0.33 
+0.34 


total models 




-0.25 
+0.42 


+0.23 
-0.31 


+0.14 
-0.23 


-0.40 
+0.44 



Table 3: Summary of systematic uncertainties in the analysis of single and di-lepton events. Ranges given in % correspond to relative 
variations around the central value. 
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amount suggested in ||14|. Varying the B hadron composition was found to produce 
negligible effect. 

the modelling uncertainty related to different semileptonic decay models has been 



calculated according to |T^. Thus the ISGW and ISGW** models have been used 
as conventional references for evaluating the semileptonic decay model uncertainty 
and this variation represents the dominant source of systematic uncertainty. 
• the finite statistics used in the simulation was checked to introduce a negligible 
systematic error. 

The summary of systematic uncertainties is given in Table ^. 

In conclusion from a fit to single and di-lepton events from data collected with the 
DELPHI detector in 1994 and 1995, the semileptonic branching fractions BR(6— 
BR(b— ^c^i'"'"), BR(b^c^£~) and the average b mixing parameter x have been measured: 



BR(6^r) = (10.71 ± O.ll(stat) ± 0.26{syst)-°oll{model))% 

BR(b^c^£+) = {8.05±0.39{stat)±0.38{syst)tofi{model))% 

BR(b^c^r) = {lM±0.35{stat)±0.25{syst)to,lt{model))% 

X = 0.126 ± 0.013(stat) ± 0m6{syst) ± 0m4{model) 



7 Analysis II: Measurement of semileptonic b decays 
from single leptons and jet-charge 

In this analysis a sample of b enriched events was obtained by applying 6-flavour tagging 
separately to each hemisphere of the event, only events with the thrust axis contained 
in the region \cos9thrust\ < 0.8 were used. The b tagging algorithm exploited only the 
information from the impact parameters of the tracks from charged particles assigned to 
the hemisphere: the cut selected 69.2 % of 66 , 12.9 % of cc and 1.1 % of uds events, so 
that the fraction of b events in the sample was Vb = 84.0%. Leptons were selected from 
all the charged particles with momentum p > 3 GeV/c, lying in the hemisphere opposite 
to the 6-tagged hemisphere within the acceptance of the HPC or muon chambers. 

The lepton was then used as a seed to reconstruct the position of the b decay vertex, 
by applying the algorithm originally developed for lifetime and oscillation measurements 
(for details, see e.g. A vertex was found in 92.5 ± 0.2 (92.3 ± 0.1)% of the cases in 

the data (simulation). The direction of the 6-hadron was then obtained by averaging the 
direction of the jet containing the lepton with the one of the vector joining the primary 
to the secondary vertex: when the vertex was not reconstructed, only the jet direction 
was used. The energy of the b hadron was computed from the sum of the energy of the 
charged and neutral particles assigned to its jet and the missing energy in the hemisphere 
(computed as described in ||21|). The resolution was a{EB)/EB — 12%. This allowed the 
entire 6-hadron four-momentum to be reconstructed, by assuming an average mass of ~ 
5.3 GeV/c2. 

Leptons from direct decays were then separated from the other sources of lep- 
tons by means of kinematics and charge correlation, as described in the following. The 
momentum of the lepton in the 6-hadron rest frame, k*, was computed by boosting back 
the lepton into the 6-hadron rest frame: the resolution was about 0"^. ~ 200 MeV/c. The 
k* spectra for b-^i^ , b-^c^i^, c— decays in the simulation were tuned as described 
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in section and varied according to the prescriptions already described to compute the 
systematic uncertainty. 

The charge of the lepton, Qe, was compared to the one of the b jet measured in the 
opposite hemisphere, Qb. Neglecting mixing, the product Xq = Qe ■ Qb should be, in 
case of perfect measurement, -1/3 (+1/3) for leptons from direct (cascade) decays. The 
charge of the b quark was determined in each hemisphere by properly combining several 
quantities (jet charge, vertex charge, charge of any kaon or lepton from b decay, charge 
of leading fragmentation particles: a detailed description of the method can be found in 
||2^ ), such that Ag actually ranged between -1 (mostly 6— >£~) and +1 (mostly 6— *c— >£^). 
Figure || shows the Xq distribution for the data and simulation. The fraction of wrong 
charge assignment, for a given Xq range, depends on several quantities related both to 
the b hadron production and decay mechanisms {B mixing, fragmentation, lepton and 
K production in b decays, b charged multiplicity, etc.) and to the detector performance 
(tracking, vertexing, particle identification), which are in some cases not well known. To 
reduce the systematic uncertainty, the fraction of correct tags was determined in the data, 
as explained in Section [7T| . 

For the previous analysis the charge correlation was only available for the di-lepton 
sample whereas Xq can be determined for all events: it should be noted however that 
the discrimination power of this variable is smaller. Therefore the two analyses are 
complementary. Only 1994 data were used for this analysis. 




Figure 2: Distribution of the charge correlation variable Xq = Qe ■ Qb for data and 
simulation. 
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7.1 Determination of the branching fractions 

The b semileptonic branching fractions were obtained by means of a binned fit- 
Leptons in the data and in the simulation were collected in two-dimensional bins, accord- 
ing to their k* and Xq values, so as to exploit fully the discriminating power of the two 
variables. The k* bins had adjustable widths, defined such as to correspond to at least 
40 entries in each bin. The range of the Xq values was divided into an even number 
(A^Aq) of bins of the equal width, 4 Xq and 25 k* bins were used. 

Events in the simulation were assigned to one of the seven classes described in Section 
|6TT| depending on their origin. Leptons from classes (d) to (g) were normalized to the 
data according to the number of hadronic events, known branching ratios and efficiency 
correction factors. The normalization factors for the classes (a), (b) and (c) were instead 
determined from the fit and used to compute the branching fractions for the direct (6— 
and cascade (6— i^c— b^c^i~) semileptonic decays. Figure | shows the fitted k* 
distribution in four different Xq bins. 

The fraction of correct charge tags in each Xq bin was determined while performing 
the fit. For this purpose, the total number of simulated events belonging to the class a 
(a=a,b,c) and falling in the i^^ (j*'*) k* {Xq ) bin {■N'^jQii, j)) were multiplied by a linear 
correction factor: 

= Mmc^{i,j) ■ (1 + 5,") 

where Af°'{i,j) is the number of data events in the same bin. The 6 coefficients would be 
zero if the simulation described the data perfectly. They were left as free parameters in 
the fit with the following constraints: 

• for a given Xq bin, 6 does not depend on k* 

• Sj = 6j = 6}^, where k is the Xq bin with opposite charge with respect to j {k = 

• = E„ -^Mc(^'i) for every a 

The first requirement follows from the fact that the Xq value is computed in the hemi- 
sphere opposite to the lepton, and is therefore uncorrelated with the value of k* and with 
all other lepton decay properties. The second constraint expresses the fact that leptons 
from direct and cascade decays populate mainly cells that are symmetric with respect to 
Xq . The third constraint ensures that the total number of events is conserved. Values of 
6 of about -7% and +4% have been obtained for classes (a) and (b,c), respectively. The 
fit results did not change significantly if the same correction was applied to the simulated 
leptons of the other classes (d-g). 

The procedure was performed separately for muons and electrons: consistent results 
were found. The psr degree of freedom was 0.95 for muons and 1.23 for electrons. 
There was no appreciable difference in the when using different models to describe the 
lepton spectra. 





BR(6^r) BR(b-^c-^£+) BR(b^c^r) 


BR(6-^r) 

BR(b-^c->£+) 

BR(b^c^r) 


1.00 0.017 -0.228 
1.00 -0.928 
1.00 



Table 4: Correlation matrix of statistical uncertainties in Analysis II. 



15 





• • • data 




simulation 


DELPHI 


[/Wl b-^l 




b— c— 1 








0.5 1 1.5 2 2.5 3 

k* GeV/c 



o 

O 

• *— ( 
U 




0.5 



1 1.5 2 2.5 3 

k* GeV/c 



^ 800 
I 700 
600 
1 500 
I 400 
I 300 
200 
100 








0.5 1 1.5 2 2.5 3 

k* GeV/c 



^ 600 

O 500 

CM 

400 

CO 

1 300 
^ 200 



100 








0.5 1 1.5 2 2.5 3 

k* GeV/c 



Figure 3: Distribution of the lepton momentum in the b rest frame, k* , shown in different 
Xq bins. 



Error Source 


Range 


ABR{b^i-) 

10-2 


ABR(b^c^£+) 
10-2 


ABR(b^c^£-) 

10-2 


electron efficiency 


±3.% 


^0.15 


^0.12 


TO.09 


misidentified electrons 










and converted photons 


±8.%, ±10% 


±0.01 


TO.03 


TO.08 


II efficiency 


±2.5% 


^0.17 


TO.09 


TO.07 


misidentified /i 


±6.5% 


<0.01 


< 0.01 


TO.07 




±9% 


±0.14 


±0.10 


±0.03 


^uds 


±22% 


±0.03 


±0.02 


<0.01 


£-btag correlation 


±1.% 


^0.05 


T0.11 


TO.03 


Rb 


0.21643 ±0.00073 [17] 


<0.01 


<0.01 


<0.01 


Rc 


0.1694 ± 0.0038 [17] 


±0.01 


±0.01 


=fo.oi 


binning 


± 2 bins 


± 0.05 


±0.05 


±0.05 


total experimental 




±0.28 


±0.22 


±0.16 


XE{b) 


0.702 ± 0.008 [14] 


< 0.01 


< 0.01 


< 0.01 


Xe[C) 

BR{b ^T^i) 


0.484 ± 0.008 [14] 
(l-28lait) [14] 
(0.459 ± 0.071)%, [2] 


=F0.02 
±0.03 
^0.01 


±0.02 

+0.20 
-0.11 

TO.04 


< 0.01 

—0.23 
+0.13 

TO. 10 


T~l I ^ / 7 Til f)-V- f) — \ 

Bll(o — ^ J/%p — >• t^t ) 

T3T3 ( r. 0+\ 

Dtxyc — >i^) 
g cc 


(0.07 ± 0.01)% [2J 

_L n Qo^o/ [1 vi 
[y.oo ± \J.oZ) /c [L(\ 

(3.19 ± 0.46)% [17] 


=F0.02 

<0.01 


±0.01 

< U.Ui 
<0.01 


T0.02 

-i-n no 
TU.Uz 

<0.01 


g ^ bb 


(0.251 ± 0.063)% [17] 


<0.01 


<0.01 


<0.01 


total systematics 




±0.28 


±0.28 


±0.27 


Semilept.mod.b — > -^[14] 
Semilept.mod.c — > ^[14] 


ACCMM (^S„) 


-0.33 
+0.53 
-0.08 
+0.06 


-0.27 
+0.44 
-0.22 
+0.09 


+0.56 
-0.84 
+0.07 
-0.05 


total models 




-U.34 
+0.53 


-U.35 
+0.50 


+U.56 
-0.84 



Table 5: Summary of systematic uncertainties in the analysis of lepton vs jet charge. Ranges given in % correspond to relative variations 
around the central value. 



17 



The final results, averaged between electrons and muons, are: 

BR{b^r) = {10.78 ±0.U{stat)±0.28{syst)-°Qll{model))% 
BR{h^c^t) = ( 7.59 ± 0.69{stat) ± 0.28(s2/st) J:^^(mo(ieO)% 
BR{h^c^r) = ( 2.00 ± 0A9{stat) ± 0.27 {syst)toil{model)% 

The average correlation matrix for the statistical uncertainties is shown in Table |[ The 
breakdown of the systematic uncertainties for the fit is presented in Table |^. The variation 
of the k* resolution causes small differences in the bins population which are included in 
the binning error. 



8 Analysis III: Measurement of semileptonic b decays 
by applying a multitag method 

A measurement of BR(6 fi) and BR(6 c(c) /i) using data collected with the 
DELPHI detector between 1992 and 1995 is presented here. Muons were identified as 



described in Section 5.1 



In this analysis the contributions of uds, c and b fiavours were separated in an inclusive 
way using a multitag method which used almost all the hadronic events, because it was 
based on a fiavour deconvolution without the need for any further cuts. One important 
by-product of the method was a systematic and independent analysis of the muon back- 
ground; as this study cannot be simply applied at electrons due to the presence of photon 
conversions, all the analysis has been performed with muons only. 

The selection of the hadronic events was the same as in Section || except that five 
charged particles instead of seven were required to select the event, and the event thrust 
axis was required to satisfy | cos^f^l < 0.75. 

The total numbers of selected events both in real and simulated data are shown in 
Table |. 

8.1 Flavour tagging 

The uds, c and b events were separated using the multivariate analysis which was pre- 



viously applied to the Ti,i/Thad determination IT^]. In each event hemisphere defined with 
respect to the thrust axis, a set of discriminating variables, called discriminators, were 
calculated, using lifetime information and event shape variables. These were combined in 



the multivariate flavour tagging algorithm [|^] and the flavour confidence algorithm ||TT 



The outputs of these two algorithms were then combined as in [jl^]. By applying cuts 
to the combined discriminator and, as in ||1 1|| , using the enhanced impact parameter tag 
to define the b-tight category, each hemisphere was classified in one of the following six 
categories: uds-loose, uds-tight, charm, b-loose, b-standard and b-tight, numbered from 
1 to 6 respectively. 

The 6 hemisphere categories provide 21 corresponding event categories and hence 21 
equations from which the 18-3 independent probabilities, ej, of classifying a hemisphere 
of flavour j in category i (j = b,c,uds and i = 1,...,6) and the 3-1 independent Rj 
values, the fractions of flavour j hemispheres in the whole sample, might be determined 
from a flt to the data. But in practice, because of a rotational ambiguity in the system, 
3 additional inputs have to be given. As in |[ll|, these were chosen to be Rc and the 
probabilities £h-tight ^l-ught of classifying charm and uds hemispheres in the b-tight 
category. 
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1992 


1993 


1994 


1995 


Total 


Simulated 
Real data 


1369 156 
486 357 


1232 678 
471 437 


2 275 552 
971 448 


712 868 
467809 


5 590 254 
2 397051 



Table 6: Total numbers of selected events for Analysis III 



1 



In this analysis the main output of this step is the determination of the probabilities 
and hence the flavour content of the different hemisphere categories, rather than 
that of Rh- The cuts on the combined discriminators have therefore been re-optimized 
with respect to ||lT|. The cut on the extended impact parameter tag, however, was 

s 

tight 



kept unchanged in order to keep the values of sl-uqht sl-tiqht unchanged from those 



'b— tight 

determined in ||11|. The value of Rb obtained was Rb = 0.21741 ± 0.00065 (stat). 

The two main features of this method are the minimal correlation between hemispheres 
(because the event vertex was computed independently in each hemisphere) and the direct 
measurement of the tagging efficiencies and of the flavour composition from the data. 
Since 1994, due to the introduction of double sided silicon detectors measuring z as well 
as r0, a better 6-flavour tagging has been achieved. 



8.2 Flavour deconvolution 

The aim of the flavour deconvolution was to extract the spectra of the muon variables 
p, pI"' and for each flavour, where p is the momentum of the muon candidate and 
pI"' and are its transverse momentum with respect to the jet axis including (pj") or 
excluding the muon in the deflnition of the jet. Hereafter any of these variables 

will be referred to as z. The inputs to the flavour deconvolution were the distributions 
of these variables for each of the six categories defined in the previous section: the 
category assigned to an identified muon was the category found by the tagging in the 
opposite hemisphere, in order to avoid correlations between the hemisphere tagging and 
the presence of the muon. 

A was then constructed using the number n^{z) of identified muons in a given 
category, i, in an interval of z: 

where Nhem is the total number of hemispheres, Rj and are the flavour fractions and 
tagging probabilities extracted from the data as just explained above, and Dj{z) is the 
spectrum of the z variable for flavour j extracted from the flavour deconvolution. The 
above formula neglects correlations between the hemisphere tagging and muon selection 
efficiencies in opposite hemispheres. 

The minimization of this function leads to a set of three linear equations for each 
z bin, where the three unknowns are the components of the spectrum in each fiavour: 
-^nds(^)' ^ci^)j ^bi^)- These quantities, and their errors, were computed by solving 
these equations. 

Thus, as a result of the deconvolution, the spectra of identified muons in the different 
fiavours were obtained. They can be written as a function of the different sources of 
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muons: 

n^ci^) = Nf,,^R,D^{z) = (^) + n^s'^(^) (2) 
n^,{z) = NhemRbD^iz) = nl^{z) + n'^z) 

where n''udsi^)^ n^J''^{z) and nj^^{z) are the distributions of background muons for different 
flavours, and n^^{z) and n\^{z) are the distributions of prompt muons coming from c and 
h decays respectively. 

This method of flavour deconvolution can also be applied to other kinds of particles 
and observables. For example, the deconvolution can be applied to all charged particles. 
The distributions obtained with charged particles are interesting results in themselves, 
but are here used only to compute the backgrounds n^^'^{z) and nj^^{z) from ri'ufgiz), as 
described in the next section. 



8.3 Background extraction and hadron misidentification proba- 
bility 

In this analysis, a background muon was deflned as any particle identifled as a muon 
that either was not a muon, or was a muon but from a light hadron (mainly pion or 
kaon) decay. Following this definition, all identified muons in uds events were taken as 
background. The misidentification probability, rjuds, was then defined as the fraction of 
charged particles identified as muons in uds events: 

Vudsiz) = ^ (3) 

where n^^dsi^) the spectrum of charged particles with the same kinematic cuts as the 
muons in the uds sector. 
This can be expressed as: 

Vudsiz) = v''{z)f:A^) + v''{z)fL{^) + V'{z)f!:dsi^) + V''iz)f^,,{z) (4) 

where ri'"{z) and rj^{z) are the misidentification probabilities for pions and kaons, f^^si^) 
and fudsi^) the fractions of pions and kaons for the uds fiavour, f^dsi^) the fraction 
of muons coming from tt and K decays in fiight and r]^{z) is their identification efficiency, 
and fudsi^) ri"{z) are respectively the fraction and the misidentification probability 
of other charged particles, which are mainly protons. The fractions for the different 



fiavours and particles have been measured in DELPHI and agree with the predictions 
obtained with the JETSET simulation program and used in this analysis. The specific 
misidentification probabilities (rj'^^z), ri^{z), ...) were supposed to be fiavour independent 
but, since the fractions of these particles are not the same in uds, c and b events, a different 
misidentification probability was evaluated for each fiavour {rjuds, Vc and rjb). Equation 
(I) was used to extract ri'^{z), taking rjuds{,z) from the data and ax-K = (z) / rj'^ (z) , rj'^'^z) 
and r]°{z) from the simulation. Then, from equations analogous to written for c and 
b fiavours, rjc and rjb were calculated. 

The misidentification probabilities obtained with this method were compared with 



those obtained using a tight anti-6 cut in Section |5TT| , and good agreement was observed. 

Once the misidentification probability for each fiavour was computed, the numbers 
of background muons per hemisphere for a variable z, i.e. the n'^^^i^z) in @, were ob- 
tained by multiplying them by the number of charged particles per hemisphere for each 
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flavour. Subtracting these contaminations from the muon candidates per hemisphere, it 
was possible to determine the distributions of prompt muons. 

8.4 Fitting of prompt muon distribution 

In order to measure the branching fractions BR{b — > /i) and BR{b c(c) — > fi), the 
following function was then minimized: 

- (nr(z.)-nr'^(z,))' 

X =J2- pmt^ — - (5) 



i=l 



where m is the number of bins, nl'^(zi) is the distribution of prompt muons measured as 
described above, and nl^' (zi) is a model expectation which can be written as: 



n 



pfi,th 
b 



(z) = NhemRb{l + BR{g^bb)) 

X [eb-,^{z)Pb^^{z)BR{b fi)+ eb^c{c)~*fi {z)Pb-,cic)~.f.{z)BR{b ^ c(c) ^ /i)](6) 

where BR{b fi) and BR{b c(c) fi) are the only unknowns, and Pb^^{z) and 
Pb^c(c)^^i{z) are the true spectra of muons coming from b fi and b c(c) /i 
decays which were taken from different models: for the central value, the ACCMM model 
has been used for 6 — > yU decays and the ACCMMl model for c — > /i decays. The 
additional terms n'^^^^^{z), n^_^j^_^_^^{z) and Ug^^^^^^z) are the contributions to prompt 
muons coming from b ^ t ^ fi, b ^ J/ip ^ fi and g ^ cc ^ ^ decays, respectively. 
The shapes of these distributions have been taken directly from the simulation, but the 



recommendations of [|1J] have been followed for their normalizations. 

The factors e^^^ and eb^c(c)^/i are global efficiency factors which contain the product 
of the efficiencies for the momentum cut (p > 3 GeV/c) and the muon geometrical 
acceptance, evaluated for each of the two considered channels, and the muon identification 
efficiency. 



8.5 Results and systematic errors 

The semileptonic branching fractions were obtained minimizing the binned of equa- 
tion (H). In order to check the validity of the method, a test was performed using simulated 
data. Figure ^ shows a comparison between the muon distributions at generation 
level and after deconvolution. A small discrepancy is visible in the b sample. The dif- 
ference between the generated values of the semileptonic branching fractions and the fit 
results were found to be 0.8% and 1.4% for the direct and cascade muons, respectively. 
These differences take into account the approximations used in the analysis. They were 
used to correct the results obtained with data and were also taken as systematic error 
contributions. 

The results obtained applying the fitting procedure to the real data are shown in 
Table ^ It can be seen that some variables, which separate the different contributions 
in different regions, are more discriminant than others. For the transverse momentum, 
b — > c(c) — > fi events are concentrated at low values, while b ^ jj, events are mainly 
situated at high transverse momentum. On the other hand in the p distribution, in the low 
momentum region both contributions are of similar importance. Thus the errors on the 
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Figure 4: Separation of spectra of candidate muons between the three flavors. The 
upper part of each plot compares the results of the deconvolution in simulated data 
(points) with the generated spectra (solid line); the lower part shows the ratio between 
these two distributions. 
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b ^ jj, 
(%) 


— > c{c) — > /i 

(%) 


X /dof 




p 


10.78 ± 0.28 


9.22 ± 0.46 


25.38/27 


1992 




10.79 ± 0.25 


9.68 ± 0.42 


25.20/32 




Pt 


iU.Yo ± U.zz 


9.81 ± 0.37 


OO /OO 

22.75/32 




P 


10.77 ± 0.29 


9.24 ± 0.50 


30.62/27 


1993 


pT 


10.68 ± 0.25 


9.77 ± 0.45 


30.02/32 




^out 

Pt 


1 n ciQ _i_ n OO 


y. lo ± U.4U 


/II ^^O /OO 

4i.DZ/oz 




p 


10.77 ± 0.18 


9.60 ± 0.25 


43.05/27 


1994 


pf 


10.73 ± 0.16 


9.43 ± 0.28 


27.74/32 




(y.OUt 
Pt 


10.62 ± 0.14 


9.54 lb 0.24 


37.16/32 




P 


10.76 ± 0.29 


9.69 ± 0.45 


18.82/27 


1995 


pin 


10.72 ± 0.24 


9.86 ± 0.41 


24.21/32 




Pt 


10.67 ± 0.21 


9.93 ± 0.36 


39.26/32 



Table 7: Fit result for the real data (the errors are only statistical). 



semileptonic branching fractions extracted using the transverse momentum distributions 
are expected to be lower than those obtained using the momentum distribution. 

Once the b semileptonic branching fractions have been fitted, it is possible to calcu- 
late the b ^ jjL and the b — > c(c) jj, spectra using the model spectra Pb^^{z) and 
Pb^c{c)^^i{z). The results are displayed in Figure |^ for each year of data taking. The 
small contributions coming from the 6 r — > /i and 6 — > J/^ — > /i decay channels, taken 
directly from the simulation, are also shown. 

Sources of systematic uncertainties have been grouped into several different categories. 
Here we comment briefly on the features that are specific to this analysis: 

• muon misidentification: The independent determination of the background distri- 
butions in this analysis is affected by 

- the values of f^, and which are the fractions of pions, kaons, muons 
(coming from vr and K decays in flight), and other charged particles in b events; 
the central values were taken from JETSET and the errors (a) in the table are 



taken from |^]; 2a ranges are taken to conservatively cover the degree to which 
the DELPHI data HI corroborated the JETSET values. 
- the misidentification probabilities specific to the particles such as rf ^ which has 
been evaluated from 'quds, the ratio ax-w, which has been taken from simula- 
tion, and rj'^ and 77", whose contribution is small and has also been taken from 
simulation. 

hemisphere tagging: in order to use the multivariate method, three parameters had 
to be fixed externally: Rc and the probabilities Sb'lught ^md £b-tight'i the variations 
of the latter probabilities correspond to their systematic uncertainties as evaluated 
in ||Tl|. The variation corresponding to the difference between the Rb value resulting 
from this analysis and the reference value used from the other three analyses was 
found to be negligible. 

analysis method: here the effects of different choices made in our analysis are con- 
sidered, namely (i) the choice of the variable (i.e. p, pj" or p°"*), (ii) the effect of 
using a looser muon selection, (iii) the infiuence of changing the number of bins of 
our variables, and (iv) the effect of the bias shown in Figure || and discussed above. 
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Figure 5: Comparison of the distributions of prompt muons for the b flavour in real 
data (dots) with the distributions obtained using the semileptonic branching fractions 
(histograms). The contributions of different processes are displayed. 
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For each year the results obtained with the three variables were averaged assuming 
complete correlation in the statistical error. After averaging over the four years, taking 
into account the correlations between the systematic errors, the results are: 

BR{b^fi) = (10.71 ±0.11(stat) ±0.28(s?/st);°i^(mode/))% 
BR{b c(c) ^ /i) = ( 9.62 ± 0.19(stot) ± OAl{syst)toAUmodel))% 

9 Analysis IV: Measurement of semileptonic b decays 
from inclusive 6-hadron reconstruction and charge 
correlation 

In this analysis the charge correlation between the b quark and the lepton produced 
in its decay was used to measure the semileptonic decay rates of 6-hadrons. The two 
different cases leading to the like charges, direct decay {b-^i~) and "upper decay vertex" 
( b^c^i~), were separated on the basis of different lepton momentum regions. 

To use the charge correlation method, 6-hadrons containing a 6-quark, Hh, needed to 
be separated from those containing a 6-quark, H^. This separation was accomplished in 
four steps: 1) by isolating bb events, 2) by reconstructing the 6-hadron decay vertex, 3) by 
identifying the tracks from the 6-hadron vertex and finally 4) by estimating the hadron 
charge. The details of these four steps are described below in section |9T|.l to |9.1|.4. After 
the separation, the sign of the charge of the 6-quark and that of the lepton were compared, 
and each lepton was classified into "like-sign" or "opposite-sign" categories. The fit of 
the like-sign spectrum was performed assuming the sample was composed of b^i'and 
6— ^c— i>£~decays, whereas the opposite-sign spectrum assumed only 6— >c—>£+ decays. 

9.1 B reconstruction and separation between Hb and 
9.1.1 Event selection 

Hadronic events were selected in the same manner as described in Section ^ and the 
event thrust axis was required to be within the region | cos 9thrust\ < 0.75 to ensure a good 
6-tagging efficiency. In addition, good detector operating conditions were required for all 
detectors, including the RICH detector, used for hadron identification. Such requirements 
led to the selection of 644 792 and 223 082 events in 1994 and 1995 data taking periods, 
respectively. Each event was then divided into two hemispheres with respect to the thrust 
axis, and the combined 6-tagging algorithm described in Section || was applied to select 
hemispheres enriched in 6-hadron content. The number of tagged hemispheres which 
contain a b quark was estimated using the same technique as in Section |6T|. A slightly 
different cut on the combined 6-tagging variable was used in this analysis, obtaining 
in simulation the following c and uds efficiencies: Eh = (42.50 ± 0.06(stat))%, Ec = 
(3.01±0.02(stat))%, Euds = (0.329 ±0.003(stat))%. This led to the purity of all 6-tagged 
hemisphere being (92.6 ± 0.3{stat))%. 

For each 6-tagged hemisphere, lepton candidates were selected in the opposite hemi- 
sphere using the same criteria as in Section ^ This method avoids introducing a bias 
on the relative fraction of the different 6-hadron species in the hemispheres where lepton 
candidates were selected. 
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9.1.2 Reconstruction of the 6-hadron vertex 

In reconstructing the 6-hadron decay vertex, the rapidity method presented in reference 
25|| was used. The reference axis for the rapidity calculation was defined by the jet 



direction obtained using the LUCLUS algorithm with the transverse momentum as the 
distance between jets and the parameter djom set to 5 GeV/c. The rapidity of each charged 
and neutral particle with respect to the reference axis was calculated, the particles outside 
the central rapidity window of ±1.5 were selected as 6— hadron decay products and used 
to reconstruct the secondary vertex. A raw 6-hadron mass and energy were computed 
from the sum of the momentum vectors of the selected particles in the jet. These values 
were corrected depending on the reconstructed mass and hemisphere energy. This led 
to a relative energy resolution of about 7% for 75% of the b hadrons which constitute a 
Gaussian distribution, with the remainder making a tail at higher energies. 



9.1.3 Identifying tracks from the 6-hadron decay vertex 

For each charged particle a probability, Pi, that the particle originated from a 6-hadron 
decay rather than from fragmentation was calculated using a neural network. It took into 
account the particle rapidity and momentum, its probability to originate from the primary 
vertex, its probability to originate from the fitted secondary vertex, the flight distance 
and the energy of the hemisphere. Figure §(a) shows the comparison between the real 
data and the simulation. 



9.1.4 Classification of and Hi 

For each hemisphere, the vertex charge Qb = J^Qi^i its uncertainty aqg = 
Pi{l — Pi) were calculated by using the probability. Pi, and the charge, Qi, of each 
particle. These values, combined with the charge of the identified kaon from 6-hadron 
decay, the jet charge and the charge of the leading fragmentation particle were fed into 
a neural network to classify a 6-hadron into Hb or Hi. The jet charge was defined as: 

Qiet = ^St^P=4^, where t is the direction of the thrust axis and 'pt is the momentum 
of the track. Using simulation, the weighting exponent k was tuned to optimize the 
probability of correctly assigning the charge of 6-hadron and was chosen to be 0.6. Figure 
H(b) shows the comparison between the real data and the simulation. 



9.2 Measurements 
9.2.1 Lepton selection 

The lepton identification was performed as in Section |5|. In addition, the lepton candi- 
date was required to originate from the 6-hadron decay vertex by requiring its probability 
Pi to be larger than 0.5. 

For each selected lepton, its momentum k*, in the 6-hadron rest frame, was calculated 



using the 6-hadron four-momentum calculated in Section |9.1.2| . Since the average resolu- 
tion on /c* is 0.1 GeV/c, the k* distribution was chosen with a bin width of 0.2 GeV/c to 
reduce migration effects. 



9.2.2 Fitting and results 

The k* distributions of leptons classified as "like-sign" and "opposite-sign" were com- 
pared to the expected spectra from simulation and the branching fractions were extracted 
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Figure 6: (a) Distribution of the track probability of real data compared to the simulation 
shown in log scale: Solid (hatched) area represents the tracks from fragmentation (from 
6-hadron decay), (b) 6-hadron charge for real data compared to the simulation: dotted 
(dashed) curve represents the Hi^ {Hi). 
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by means of a binned fit. The background contributions which may arise from non- 
b events, non-6-decay products and wrongly identified leptons were estimated from the 
simulation and subtracted. Any incorrectly determined charge of the 6-quark led to the 
misclassification of leptons from like-sign to opposite-sign and vice versa. The amount 
of misclassified leptons was first estimated from the simulation and used in the fit of the 
lepton spectra. The fraction of each type of decay obtained from the fit was then used 
to adjust the amount of misclassified leptons. This process was repeated until the fitting 
results converged. 

The following results have been obtained, and Figure ^ shows the results of the fit 
using the ACCMM model, where the uncertainties are only statistical: 





1994 


1995 


combined 


BR(6^r)(%) 
BR{b ^ c ^ £+){%) 
BR(6-^c^r)(%) 


10.78 ±0.18 
8.02 ±0.31 
1.33 ±0.32 


10.67 ±0.30 
7.92 ±0.52 
1.36 ±0.50 


10.75 ±0.15 
7.99 ±0.27 
1.34 ±0.30 



The following correlation matrix was found: 





BR(6-^r) BR(b-^c-^£+) BR(b-^c^r) 


BR(6^r) 

BR(b^c-^£+) 

BR(b^c^r) 


1.00 -0.077 -0.350 
1.00 -0.603 
1.00 



9.3 Systematic uncertainties 

Since the b reconstruction and the charge evaluation of the 6-hadron were done in the 
hemisphere where the lepton candidate was found, the correlation between the lepton 
selection and the charge determination of the b hadrons must be studied. Although the 
lepton information was not included in the training of the neural network to obtain the 
charge of the 6-hadron, a small correlation of pu = 1.036 ± 0.005 was found, where pu 
represents the ratio of efficiencies to tag a hemisphere which contain a lepton over all 
hemispheres. This was used to reweight the Monte Carlo events, and twice the statistical 
error on pu was used to obtain the contribution to the systematic uncertainty. 

A more critical bias exists between the neural network output and the 6-hadron com- 
position. The neural network output for a hemisphere containing a charged 6-hadron 
was more likely to give the correct charge of the 6-quark than a hemisphere containing 
a neutral 6-hadron. The effect of this bias was to increase the likelihood of incorrectly 
determining the charge of the 6-quark for neutral 6-hadrons. However, artificially adjust- 
ing the Monte Carlo weight to account for this bias resulted in very little change in the 
branching fractions. A more critical approach was to compare the measured branching 
fractions with the ones obtained without the charge separation. Without the separation, 
the lepton spectrum contained the contributions from the direct decay and both modes 
of the secondary decays. The fit of the three modes was performed by alternatively fixing 
one rate of the two secondary decays modes, starting with the rate of b ^ c ^ i fixed to 
the result of the analysis, until the fit converged. The difference between the branching 
ratios obtained in this fit and the ones obtained with the charge separation was used as 
a systematic uncertainty. 

The contributions to the systematic uncertainties of the correlation studies are shown 
in the first part of Table ||. Other sources considered for systematic uncertainties are as 
follows: 
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Figure 7: Lepton momentum spectra in the 6-hadron rest frame. Plot (a) ((b)) shows 
the result of the fit with the ACCMM model to the like-sign (opposite-sign) sample . 
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Lepton selection: 

The muon and electron identification efficiencies and the background due to hadron 
misidentification were varied considering their measurement uncertainties in the 



data-simulation comparisons (see Sections |5TI|, ^.21 ) as in Analysis I. The residual 
contamination in the electron sample due to converted photons has been varied by 
± 10%. 
6-tagging 

The efficiencies to tag c and uds quarks, as well as the values of Rb and Ruds, were 
varied in the same manner as in Analysis I. The correlation between the lifetime tag 
and the lepton tag was found to be pe = 1.057±0.005 and = 1.041 ±0.005. These 
values were varied by twice their statistical uncertainties. 
Fitting 

The uncertainty due to the finite Monte Carlo statistics in the lepton spectrum 
fitting procedure was evaluated. 
6-hadron composition 

The production fraction for Af, was taken from [|] and set to (10.ll3j)%, and the 



semileptonic branching fraction was set to BR(Ab iu X)) = (7.4 ± 1.1)% ||26| . 
Models 

The mean fractional energy of c hadrons was varied according to ||14| . 
The lepton distribution from the "upper vertex" was studied by varying the contri- 
butions of Dg i~X and D^{D^) — > l^X as suggested in reference [T^ . 
The modelling uncertainty related to the branching fractions assumed for 6 — >^ r — >^ 
b ^ J/"^ ^ i and to different lepton decay models was also calculated according to 
0,10] and 0. 



The summary of the different contributions to systematic uncertainties is given in 
Table |^. In conclusion, with the method of charge correlation, the following results have 
been obtained from the data collected with the DELPHI detector in 1994 and 1995: 



BR{b r) = (10.75 ± 0.15{stat) ± 0.28{syst)-^oil{model))% 
BR{b ^c^r) = (7.99 ± 0.27{stat) ± 0.28(s|/st);°:?J(mode/))% 
BR{b^c^t) = (1.34±0.30(stat) ±0.27(st/st)l[5J^(mode/))% 



10 Combinations of results 

A comparison of the results obtained in the different analyses described in the previous 
sections is shown in Table 0. A procedure to combine them in order to produce a final 
set of physical parameters has been developed. The basic technique, named Best Linear 
Unbiased Estimator (BLUE) ||2^, determines the best estimate x of a physical parameter 
built by a linear combination of measurements Xi obtained by several experiments; the 
coefficients of the combination are built from the covariance matrix Eij of the measured 
quantities. The method may be easily applied to determine several physical parameters 
simultaneously, by replacing that matrix with the more general one Eiaj/3 where the 
indices i,j refer to the experiments ( here analyses I to IV ) and a, (5 identify the different 
physical parameters (here BR(6^£^) , BR(b^c^£+) etc.). 

In order to apply this technique, it is necessary to estimate the full error matrix E 
including the off-diagonal elements; it has been determined as the sum of a statistical part 
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and a systematic part with the latter accounting for the uncertainties on the parameters 
used by the analyses and obtained from other measurements. 

The statistical part has been built by splitting the statistical error aia of each pa- 
rameter a determined by the analysis i into two terms: the first one is computed from 
the observed number of leptons and is considered as fully correlated between different 
measurements; the other term is computed in order to keep invariant the total error and 
is assumed to be uncorrelated. 

The estimation of the correlation between the parameters of different analyses is more 
complicated, as it is necessary to account for the correlation already present inside each 
single analysis. A reasonable criterion for that is to build the covariance elements by 
multiplying the correlated parts of the two aia, described above, and by applying a 
correlation factor determined as an average of the correlation coefficients resulting from 
the different analyses. 

The described procedure can be applied only for identical data samples, while the dif- 
ferent analyses used somewhat different data samples; as a consequence the full statistics 
has been divided into non-overlapping subsamples and the described procedure has been 
applied to each one of them. To do this the statistical uncertainties on the measurements 
have been scaled by the ratio of the square root of the number of events used by the 
corresponding analysis and the square root of the number of events in the subsample 
itself. These subsamples do not contain any common event and may be assumed un- 
correlated; the total covariance matrix may then be obtained by summing the inverse of 
each covariance matrix and inverting again. 

A special care has been put in handling the results of the multivariate analysis which 
builds up the prompt muon distributions by a linear combination of distributions obtained 
in 6 categories; the overlap with the ^-tagged sample used by the other analyses has been 
conservatively assumed as corresponding to the category with the biggest purity and 
therefore the biggest weight. 

The systematic part of the error matrix has been evaluated by expressing a linear 
dependence on the external parameters of each result, and propagating the uncertainties 
on the parameters themselves; this corresponds to building up the sum of a set of error 
matrices, one for each uncertainty source, with correlation factors equal to 1 for all pairs 
of results affected by the corresponding external parameter, while the systematic errors 
relevant to only some of the results have been added as uncorrelated. The errors arising 
from the uncertainties on the decay models have not been used in the combination to 
obtain a result where the dependence on them is most explicit; as these errors give the 
biggest contribution to the total error this also protects from the instabilities described 



in the cited paper and in others dealing with this topic ||27| , [28|| . The total systematic 
covariance matrix thus obtained has then been summed to the statistical covariance 
matrix; the inverse of the sum has been used to weight the four analyses results and find 
the combined value along with the total error. 
The following results have been obtained: 

BR(6-*r) = (10.70 ±0.22)% 

BR(b^c^£+) = (7.98 ±0.30)% 

BR(b^c^r) = (1.61 ±0.26)% 

X = 0.127 ±0.014 



where the total error, excluding model effect, is quoted; the global of the fit is 1.52 for 
12-4=8 degrees of freedom. 
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The statistical contribution to the total error has been obtained by propagating the 
statistical uncertainties on the four analyses output to the combined values. The system- 
atic uncertainties breakdown on the combined values have been obtained by combining 
the error sets given for each analysis, using the same coefficients used to obtain the central 
values; this is equivalent to observing the effect of changing the combined values by la 



for each of the error source. The full table of errors is shown in Table 11; the correlation 



matrix for the statistical and total uncertainties is shown in Table |T2 . 

To investigate the effect of the main assumptions done in this combination ( estimation 
of the correlated part of the error, estimation of the correlation coefficient between differ- 
ent parameters determined in different analyses ) the procedure has been repeated after 
changing them slightly. The off-diagonal element in the error matrix has been changed 
using the most conservative assumption where a result does not add any information to 
another one having a smaller uncertainty. Different estimations of the correlation coeffi- 
cient between different parameters in different analyses have also been tried. Compatible 
results have been obtained. The combination performed using a covariance matrix built 
from the statistical errors only was also found to give very similar results. 



11 Conclusions 

Four different analyses have been used to measure the semileptonic branching frac- 
tions for primary and cascade b decays in hadronic Z decays from the data collected by 
the DELPHI experiment at LEP. Results are compatible and a global average has been 
obtained: 



BR(6^r) = (10.70 ± 0.08(stat) ±0.21(syst);[}j°(moc/e/))% 

BR(b-*c^£+) = {7.98±0.22{stat)±0.21{syst)toit,imodel))% 

BR(b->c^r) = (1.61 ± 0.20{stat) ± 0.17(s?/st)+°:^°(morfe/))% 

X = 0.127 ±0M3{stat)±0m5{sy St) ±OmA{model) 



The present result is compatible with and more precise than the previous DELPHI one 
1^. It hence supersedes it. It is also compatible with the recent results of the semileptonic 
branching fraction obtained at LEP and with theoretical calculations 
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A Appendix 

A.l Single lepton likelihood 

The first part of the likehhood was constructed assuming a Poisson probabihty, using 
the single lepton spectra in data and simulation, subdivided in 25 x 25 bins in the {pt,Pi) 
plane. The bins were chosen in such a way to have approximatively the same amount 
of data in each bins. Nine classes were used, corresponding to the classes (a) to (g) 



mentioned in section |6T| , with classes (f) and (g) splitted in two, for bb and non-66 
events. 



C, = ln{L,) = Y.Y1 {DATit,j)lniE{t,j)) - E{t,j)} 

i=l j=e,ti 
N 1 

class 

m,j)= 5^m«)MC(^,j,«)} 

«=i 

where DAT{i,j) represent the data and MC{i,j) the simulated spectra, respectively. 
The V{a){a = 1,3) coefficients are the ratio between the unknown branching fractions 
and the corresponding values used in the simulation: 



, BR(6^r) BR(b^c^£+) BR(b- 



BR{b^i-),,J ^ ^ BR(b^c^£+)_' ^ ^ BR(b^c^£-)_ 

whereas the V coefficients corresponding to lepton classes (d) to (g) are fixed to the values 
given in Table 13. 



A. 2 Di-lepton likelihood 

The second part of the likelihood was constructed assuming a Poissonian probability, 
using the di-lepton spectra in data and simulation, subdivided in 7x 7 bins in the combined 
momentum variables (pj!**"^,/)^"^'). 

The bins were chosen in such a way to have approximatively the same amount of data 
in each bins. Twenty classes were used, according to the different possible combinations in 
the two opposite hemispheres of the single-lepton classes (a) to (g) mentioned in section 
6T|. 

C2=ln{L2) = ^ ^ { DATsame{i,j)ln{Esame{i,j))-Esame{i,j) + 
i=l j=ee,fifj,,efi 

DAToppXi, j)ln{EoppXi, j)) - EoppXi,j)} 

-^'^ class 

Esaraeihi) = ^ {S (a) MCsame{i , j , a)} 
AT 

^^'^ class 

Eopp{i,j)= Y,{0{a)MCopp{i,],a)} 

a=l 

where DATgameih j) ( DAToppXh j) ) represent the spectra of di-leptons in data, in op- 
posite hemispheres, having the same (opposite) charge and MCgameih j) (MCoppXh j)) 
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represent the simulated spectra. The S{a) ( 0{a) ) coefficients depend on the ratio 
between the unknown branching fractions and the corresponding values used in the sim- 
ulation and on the mixing probability x- For example for the first and the second classes, 
containing {b^£~ ,b^i~) and {b^i~ ,b^c^i^) di-leptons, respectively: 



5(1) = 2x{l-x)V{iy = 2x{l-x){ 



■BR{b^i-)s 

0(1) = (l-2x(l-x))nif = (l-2x(l-x))( 



2 n n.-.n BR(6^r) 



5(2) = (l-2x(l-x))P(l)P(2) = (l-2x(l-x)) 
0(2) = 2x(l-x)P(l)^(2)=2x(l-x) 



BR(6^^ )sim 

BR(6^£-)BR(b^c->£+) 



BR(6^£-),,^BR(b^c^£+), 
BR(6^r)BR(b^c^£+) 



BR(6^£-),,„BR(b^c^£+)_ 
The total likelihood is the sum of the single and the di-lepton likelihoods: 

£ = £i + £2 

In the fit -P(l), P(2),P(3) and x are free parameters, whereas the V coefficients corre- 
sponding to lepton classes (d) to (g) are fixed to the values given in Table 
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bource 


A(6 ^ /i) 


A{b — > c(c) /i) 


muon efficiency (±2.5%) 


T 0.190 


T 0.182 


K (±2a) 
A" (±2^) 

n (±2^) 


T 0.004 
T 0.002 
± 0.003 
T 0.001 


T 0.008 
T 0.007 
± 0.009 
T 0.001 


r]^ (±2(t) 
ai^^ (±2o-) 
r/^ (±2(t) 
r]" (±2a) 


T 0.022 
± 0.008 
T 0.004 
T 0.001 


T 0.120 
T 0.035 
T 0.004 
T 0.001 


= 0.2170 ±0.0009 
Rc = 0.1734 ±0.0048 

ef-%kt (±15%) 
et_,,^,, (±7%) 


< 0.01 

< 0.01 
± 0.023 
± 0.007 


< 0.01 

< 0.01 
± 0.010 
± 0.028 


Variable 
Muon quality 
Binning 
Bias of the method 


± 0.080 
± 0.082 
± 0.078 
± 0.080 


± 0.150 
± 0.082 
± 0.079 
± 0.136 


MC statistics 


± 0.088 


± 0.163 


XE{b) = 0.702 ±0.008 


± 0.093 


± 0.165 


BR{c ^e) = (9.85 ± 0.32)% ||l3 
BR{b ^r^£-) = (0.459 ± 0.071)% [|| 
BR{b J/ip ^ £-£+) = (0.07 ± 0.01)% HI 
5i?(^ ^ cc) = (3.19 ±0.46)% 17 
BR{g bb) = (0.251 ± 0.063)% 17 


T 0.001 
T 0.014 
T 0.018 
± 0.009 
T 0.033 


T 0.002 
T 0.096 
T 0.011 
± 0.010 
T 0.043 


total systematic 


± 0.28 


± 0.41 


b i ACCMM+ll^^^^ 


-0.35 
+0.43 
-0.11 
+0.11 


+0.52 
-0.48 
-0.12 
+0.02 


total models 


+0.44 

-0.37 


+0.52 
-0.49 



Table 8: Analysis III: Systematic uncertainties (%) for BR{b — > /i) and 
BR{b c(c) ^ fi) 
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Source 



Range 



xlO-2 



ABR 

xlO-2 



ABR 

xlO-2 



^-charge tag correlation 
NN bias on the 6-charge 
6-hadron composition 
electron efficiency 
muon efficiency 
Misidentified e 
Misidentified /i 
Converted 7 



±1% 
see text 
see text 



^0.08 
^0.08 
^0.04 



^0.03 
^0.15 
^0.02 



^0.09 
^0.11 
^0.04 



±3% 
±2.5% 

±8% 
±6.5% 
±10% 



^0.18 
^0.13 
±0.01 
±0.01 
±0.01 



^0.04 
^0.05 
^0.11 
^0.08 
^0.04 



TO. 15 
TO. 10 
TO.08 
TO.05 
TO.03 



^uds 

l-h tag correlation 
Rb 

Rr 



±9% 
±22% 
±1% 
0.21643 ±0.00073 
0.1694 ±0.0038 



17 



17 



<0.01 
<0.01 
TO.09 
<0.01 
<0.01 



TO.01 
±0.01 
TO.03 
<0.01 
<0.01 



TO.01 
TO.01 
TO.09 
<0.01 
<0.01 



MC statistics 



TO.03 



TO.01 



TO.03 



XEip) 

xe{c) 

b^W^Ds 

BR(6 ^ r ^ 
BR(6 ^ J/^ 
BR(c ^ tj 
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0.702 ±0.008 

0.484 ±0.008 
(1.28tS|?) 
(0.459 ±0.071)% H 

(0.07 ±0.01)% |] 
(9.85 ±0.32)% IP 



±0.03 
TO.01 

+0.04 
-0.04 

TO.02 
TO.02 
TO.01 



±0.05 
±0.01 

-0.09 
+0.08 

TO.07 
±0.01 
TO.05 



±0.07 
TO.01 

+0.03 
-0.03 

< 0.01 
TO.01 
TO.02 



Total systematic 



±0.28 



±0.27 



±0.28 



Decay models 



b ^ £ model 
c — ^ £ model 


ACCMM {Vj^^ZJ 
ACCMMl tfa^cM^l) 


-0.23 
+0.42 
-0.07 
+0.07 


+0.36 

-0.58 
+0.06 
-0.05 


+0.04 
-0.04 
-0.21 
+0.09 


Total Models 


-U.24 
+0.43 


+0.36 
-0.58 


-U.21 
+0.10 



Table 9: Analysis IV: Summary of systematic uncertainties. Ranges given in % corre- 
spond to relative variations around the central value. 





Analysis I 


Analysis II 


Analysis III 


Analysis IV 


BR(6^r)% 
BR(b^c^£+)% 
BR(b^c^r)% 
(BR(b^c^£+) + 
BR(b^c^r))% 


10.71 ± 0.11 ± 0.26+0:42 
8.05 ± 0.39 ± O.SStoil 
1.64 ± 0.35 ± 0.2510.23 

9.69±0.24 ± 


10.78 ± 0.14 ± 0.28+0:53 
7.59 ± 0.69 ± 0.28~o:5o 
2.00 ± 0.49 ± 0.27Io:84 

9.59±0.30 ± 0.4ll^i^ 


10.71 ± 0.11 ± 0.28+0:4! 
9.62 ± 0.19 ± OAlt^'ll 


10.75 ± 0.15 ± 0.28+0:43 
7.99 ± 0.27 ± 0.28~o:io 
1.34 ± 0.30 ± 0.27Io:58 

9.33±0.26 ± 0.52lg:^^ 



Table 10: Comparison of the results of the different analyses. The measurements are shown using boldface characters, whereas slim-face 
characters are used for sums which are only shown for comparison. The first uncertainty is statistical, the second is systematic and the 
third is due to the uncertainty on the semileptonic model. 



00 



Error Source 


Range 


ABR(6^r) 
10-2 


ABR(b^c^£+) 

10-2 


ABR(b-^c-^£-) 

10-2 


Ay 

A, 

10-2 


statistical 




^0.08 


^0.22 


TO. 20 


±1.3 


electron efficiency 




±3% 


TO.09 


TO.08 


TO.04 


±0.01 


misidentified e 




±8% 


^0.02 


^0.05 


TO.03 


±0.04 


converted photons 




±10% 


<0.01 


=F0.02 


<0.01 


=F0.03 


/i efficiency 




±2.5% 


^0.15 


^0.12 


^0.04 


TO.01 


misidentified /i 




±6.5%; 17% 


<0.01 


^0.03 


^0.03 


TO.07 


£c 




±9% 


±0.01 


±0.01 


±0.03 


±0.02 


^uds 




±22% 


±0.01 


<0.01 


<0.01 


<0.01 


i — b correlation 




±1% 


^0.03 


^0.05 


^0.02 


TO.02 


other sources 




±0.09 


±0.10 


±0.05 


±0.5 


XE{b) 




0.702 ± 0.008 [14] 


TO.01 


±0.03 


±0.02 


±0.05 


xe{c) 




0.484 ±0.008 


^0.01 


<0.01 


<0.01 


±0.04 


b^W^Ds 




(1-28+^j?) [14] 


±0.02 


±0.08 


TO. 10 


TO.05 


BR{b ^T^i) 
BR{b J/iP ^ 


n 


(0.459 ± 0.071)% [2] 
(0.07 ±0.01)% [2] 


^0.01 
^0.02 


^0.02 
^0.01 


TO.08 
<0.01 


±0.04 
TO.06 


BR(c i) 




(9.85 ±0.32)% [17] 


^0.01 


<0.01 


TO.02 


TO.01 


g cc 




(3.19 ±0.46)% [17] 


<0.01 


<0.01 


<0.01 


<0.01 


g -^bb 




(0.251 ±0.063)% [17] 


TO.01 


^0.01 


<0.01 


±0.01 


total systematic 




±0.21 


±0.21 


±0.17 


±0.5 


Semilept.mod.b — 
Semilept.mod.c 


414] 
£[14] 


ACCMM (™„) 
ACCMMl(lig™2^ 


-0.28 
+0.44 
-0.09 
+0.08 


+0.10 
-0.02 
-0.19 
+0.07 


+0.37 
-0.47 
+0.05 
-0.04 


-0.3 
+0.3 
-0.3 
+0.3 



Table 11: Systematic uncertainties associated to the combined results; the effect of sources relevant to only one analysis has been 
summarized in a single value labelled "other sources". 



CO 



□ 



40 





BR(6^£-) BR(b^c^£+) BR(b^c^£-) x 


BR(6^r) 

BR(b^c^^+) 

BR(b^c^r) 

X 


1. -0.066 -0.051 0.018 
0.545 1. -0.733 -0.091 
0.231 -0.277 1. 0.038 
0.039 -0.040 0.018 1. 



Table 12: Correlation matrix of combined results. On the upper-right side the statistical 
coefficients are reported, on the lower-left side the statistical+systematic coefficients are 
shown. 



